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ABSTRACT 

Aims. We perform an analysis of the X-ray superbubble in the N 206 H n region in the Large Magellanic Cloud using current gener- 
ation facilities to gain a better understanding of the physical processes at work in the superbubble and to improve our knowledge of 
superbubble evolution. 

Methods. We used XMM-Newton observations of the N 206 region to produce images and extract spectra of the superbubble diffuse 
emission. Morphological comparisons with Ho- images from the Magellanic Cloud Emission Line Survey were performed, and spec- 
tral analysis of the diffuse X-ray emission was carried out. We derived the physical properties of the hot gas in the superbubble based 
on the results of the spectral analysis. We also determined the total energy stored in the superbubble and compared this to the expected 
energy input from the stellar population to assess the superbubble growth rate discrepancy for N 206. 

Results. We find that the brightest region of diffuse X-ray emission is confined by a Ha shell, consistent with the superbubble model. 
In addition, faint emission extending beyond the Ho- shell was found, which we attribute to a blowout region. The spectral analysis 
of both emission regions points to a hot shocked gas as the likely origin of the emission. We determine the total energy stored in the 
bubble and the expected energy input by the stellar population. However, due to limited data on the stellar population, the input energy 
is poorly constrained and, consequently, no definitive indication of a growth rate discrepancy is seen. 

Conclusions. Using the high-sensitivity X-ray data from XMM-Newton and optical data from the Magellanic Cloud Emission Line 
Survey has allowed us to better understand the physical properties of the N 206 superbubble and address some key questions of 
superbubble evolution. 
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1. Introduction 

Massive stars, via their fast stellar winds and subsequent super- 
nova explosions, are responsible for energising and enriching 
the interstellar medium (ISM), and are the source of the ISM's 
dynamic hot-phase. Given that massive stars usually form in 
groups, their collective mechanical output into the surrounding 
ISM creates so-called 'superbubbles', 100-1000 pc diameter 
shells of swept-up interstellar material that contains a hot 
(10 6 K), shock-heated X-ray emitting gas. Superbubbles are 
one of the primary engines for driving the morphology and 
evolution of the multi-phase ISM and are of key importance in 
the understanding of matter recycling in galaxies. |Mac Low &| 
|McCray ( 1988) showed that these objects could be considered 
scaled-up versions of stellar wind blown bubbles and expanded 
on the classical analytical model of [Weaver et al.| ( fl977| ) in this 
area. More recent numerical models of superbubbles are capable 
of predicting temperature distributions and X-ray luminosities 



(see |Breitschwerdt & de Avillez 2006 Rodriguez-Gonzalez 
et al.|201 1 1 for example). 



Owing to the soft nature of superbubble X-rays (< 2 
keV) and the location of massive star forming regions in the 
Galactic disk, these objects in our Galaxy are almost completely 
obscured by foreground absorbing material. Conversely, the 
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superbubble population of the Large Magellanic Cloud (LMC) 
is ideal for study. The LMC, which is a dwarf irregular galaxy 
with indications for spiral structures, is one of the closest 
neighbours to our Galaxy. The distance to the LMC of 48 
kpc (Mac ri et al. 1 12006), its modest extinction in the line of 
sight (average Galactic foreground A^h ~ 7 x 10 20 cm" 2 ), and 
its small inclination angle of ~ 40° (Feast 1991) ) make it an 



ideal laboratory for the multiwavelength study of superbubbles. 
This allows one to trace the diffuse X-ray emitting gas and 
constrain the mass and energy input to the superbubbles using 
the observed stellar populations. 



N 206 Penize|[T956l also known as DEM L221) is an Hn 
region located in the southeast of the LMC that is excited by 
the winds of the massive stars in the young NGC 2018 cluster, 



and the LH 66 and LH 69 OB associations ( |Lucke & Hodge 
|1970[ ). [Gorjian et al.| ( |2004| > used Spitzer observations of N 206 



to identify regions of very active star formation in the shell of 
cold material surrounding the H n region. More recently, Romita 



|et al.| ( |2010| l used a multiwavelength approach to identify and 
classify many of the young stellar objects (YSOs) in N 206 and 
determined an above average star formation rate of some ~ 5 
times higher than the LMC as a whole. This Hn complex also 
harbours the supernova remnant SNR B0532-71.0 ( |Mathewson| 



& Clarke 1973b| |Williams et al.||2005 ) and contains an X-ray 



superbubble. The diffuse X-ray emission in the Hn region was 
first identified and attributed to a superbubble by |Dunne et al.| 
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Fig. 1. XMM-Newton EPIC false-colour mosaic image of SNR 
B0532-71.0 and the N 206 superbubble. Red corresponds to the 
0.3-1 keV energy range, green to 1-2 keV, and blue to 2-8 keV. 
The images have been Gaussian smoothed with a 3-pixel kernel. 



riods of high particle background, which resulted in a combined 
effective exposure of 32 ks (9 ks + 23 ks). We created mosaic 
images in the 0.3-1 keV, 1-2 keV, and 2-8 keV bands by com- 
bining the data from each MOS camera in each observation to 
produce a false-colour X-ray image of the region (shown in Fig. 
fT), which was smoothed using a Gaussian filter. The superbubble 
X-ray emission is clearly visible to the west of SNR B0532-71.0 
in the centre of the field of view (FOV). The X-ray emission is 
very soft with the vast majority of the emission being < 1 keV. 

2.2. Optical 

The Magellanic Cloud Emission Line Survey (MCELS) obser- 
vations ( |Smith et al.|2006[ ) were made with the 0.6 m University 
of Michigan/Cerro Tololo Inter- American Observatory (CTIO) 
Curtis Schmidt Telescope equipped with a SITE 2048 x 2048 
CCD, producing individual images of 1.35° X 1.35° at a scale of 
2.3" pixeF 1 . The survey mapped both the LMC (8° x 8°) and 
the SMC (3.5° x 4.5°) in narrow bands covering [Om]^5007 A, 
Ha, and [S ii]/167 16,6731 Ain addition to matched green and red 
continuum bands, which are used primarily for the subtraction of 
the stellar continuum from the narrow band images. The survey 
data were flux-calibrated and combined to produce mosaicked 
images. Cutouts from the mosaics of each LMC emission line 
image around N 206 were used for our analysis. 



(2001) using ROSAT HRI observations and Ha images from 
PDS scans of the Curtis Schmidt plates of |Kennicutt & H odge 
1986). These authors found the X-ray emission to be confined 



by an He structure, indicative of a superbubble scenario. 

In this paper we present an analysis of the N 206 superbubble 
using current generation facilities to gain a better understand- 
ing of the physical processes at work in the N 206 superbubble. 
In Sections 2 we outline the XMM-Newton and optical obser- 
vations. In Section 3 we descibe the data analysis including a 
morphological comparison to Ha emission and spectral analy- 
sis of the superbubble X-ray emission. In Section 4 we discuss 
the results of the analysis and address some key evolutionary 
questions of the N 206 superbubble. A summary of the work is 
presented in Section 5. 

2. Observations 

2.1. X-rays 

Owing to the presence of SNR B0532-71.0 in the N 206 Hn 
complex, the region has been observed with both Chandra 
( |Weisskopf et1QT|2"0"02"] l and XMM-Newton ( |Jansen et aL|2"00l] l 
with the SNR at the aimpoint. Unfortunately, the Chandra data 
are not suitable for the purposes of our study because they 
do not cover the superbubble emission region. The European 
Photon Imaging Cameras (EPIC, IStriider et al.||2001[ ITurner 



|et al.|2001| l onboard XMM-Newton have observed SNR B0532- 
71.0 on two occasions (Obs. IDs 0089210101 and 0089210901, 
PI R. Williams). Because the PN camera operated in large- 
window mode for Obs. ID 0089210101, only a fraction of the 
N 206 superbubble was covered, making these data unsuitable 
for our analysis. The PN was not in operation during Obs. ID 
0089210901. Hence, only MOS data are available for each ob- 
servation which were in imaging full-frame mode for a total of 
~ 5 1 ks during the observations. We reduced each of the datasets 
using the standard reduction tasks in S AS 10.0.0, filtering for pe- 



3. Data analysis 

3. 1 . X-ray point sources 

Source detection was performed using the SAS detection 
metatask edetect_chain across various energy bands to im- 
prove detection sensitivity. Several sources were detected in 
the region of the superbubble, which are shown in Fig. [2] with 
their multiwavelength counterparts obtained from the SIMBAD 
database given in Table [T] We detected a high-mass X-ray bi- 
nary (HXMB) XMMU J052952.9-705850 (Source 1) and a can- 
didate HMXB USNO-B1 .0 0189-00079195 (Source 3), b oth of 
which were identified by |Shtykovskiy & Gilfanov] (p005). Each 
of these sources is likely associated with the N 206 H n complex. 
Additionally, we detected an X-ray source in the NGC 2018 
cluster, namely the 04-5III(f) star HD 269676. The remaining 
sources include the foreground star HD 269669, the likely back- 
ground AGN [SHP2000] LMC 21 1 ( [Sasaki et al.|[2000l >, and a 
source without a counterpart. Detected sources located in the re- 
gions of diffuse emission were masked in the subsequent spectral 
analysis based on the optimum extraction regions determined us- 
ing the SAS task eregionanalyse. 



3.2. Morphology 

We compared the soft band 0.3-1 keV X-ray image to the 
MCELS Ha image of the superbubble to assess the superbub- 
ble morphologies, shown in Fig. [3 We found that the brightest 
X-ray emitting region of the supeibubble is confined by an Ha 
structure, as found previously by |Dunne et al. (2001 1, labelled 
as Region A in Fig. [5J;. Additionally, we identified an apparent 
blowout region to the north (labelled as Region B in Fig. |3};) 
that is surrounded by a fainter and more dispersed Ha structure. 
The presence of a blowout is also supported by the IR images 
of |Gorjian et al.| ( |2004| l, which show that the cooler material 
surrounding the Ha shell is also dispersed in the region of the 
blowout. The Ha structure confining Region A appears to be part 
of a large approximately circular Ha shell (labelled as Region 
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Table 1. X-ray point sources in the superbubble region with SIMBAD counterparts 



Source 


RA 


DEC 


SIMBAD counterpart 


Type 


1 


05:29:52.800 


-70:58:46.56 


XMMU J052952.9-705850 


HMXB 


2 


05:29:54.192 


-70:58:04.08 






3 


05:30:02.016 


-71:00:46.80 


USNO-B1.0 0189-00079195 


cand. HMXB 


4 


05:30:40.632 


-71:00:39.24 


[SHP2000] LMC211 


X-ray source 


5 


05:31:02.976 


-71:06:09.00 


HD 269669 


GO star 


6 


05:31:15.144 


-71:04:10.20 


HD 269676 


04-5III(f) star 




Fig. 2. Detected X-ray sources plotted on a cutout from the false- 
colour X-ray mosaic image in Fig. [T]in the region of the N 206 
superbubble. Red corresponds to the 0.3-1 keV energy range, 
green to 1-2 keV, and blue to 2-8 keV. See Table [T] for source 
details. 



SB in Fig. [3};) with ~ 8'diameter (corresponding to ~ 112 pc 
at the LMC distance), which encompasses the OB associations 
in the Hn region. Interestingly, no X-ray emission is detected 
from the majority of Region SB. This is somewhat unusual, as 
we would not expect the hot X-ray emitting gas to be distributed 
with such non-uniformity within the superbubble. Likely reasons 
for this are discussed in Section |4~T1 



3.3. EPIC spectra 

Spectral analysis of extended diffuse emission is complicated 
by the contributions of both the instrumental and astrophysical 
backgrounds. In the case of the astrophysical background, the 
observational FOV of our analysis is largely free of extended dif- 
fuse emission apart from the superbubble. Hence, a background 
extracted from these data should adequately account for the as- 
trophysical background, given that we can assume it to be uni- 
form over such a small area. The only problem with such a back- 
ground selection in this case is the variation in spectral response 
of the detector, which can be easily corrected for. However, this 
may not be adequate to account for the instrumental background. 
The instrumental background below 2 keV for the EPIC instru- 
ments is dominated by fluorescence lines caused by the interac- 
tion of high-energy particles with the material surrounding the 
detectors ( |Lumb et al] [20021 |Read & Ponman|[2003] ). For the 
MOS cameras these are the Al Kor and Si Ka lines at ~ 1.49 



keV and ~ 1.75 keV, respectively. Variation in the strength of 
these lines across the detector makes the background region sub- 
ject to a different instrumental background than the superbub- 
ble. One background selection method that can account for this 
effect is the use of the so-called 'blank-sky' data. The XMM- 
Newton EPIC background working group has produced these 
blank-sky data sets for each EPIC instrument over all the read- 
out modes and filters (Carter & Read 2007 ). The blank-sky data 
have been produced by merging several observations and remov- 
ing point sources to leave an observational dataset representa- 
tive of the average astrophysical background and the instrumen- 
tal background. Thus, using these blank-sky data, we can se- 
lect a background from the same region on the detector as the 
superbubble in our observational dataset, which can adequately 
account for the instrumental background and the astrophysical 
background, assuming the averaged astrophysical background 
from the blank-sky data is representative of the astrophysical 
background in our observations^ 

We assessed the effect of each of these background selec- 
tion techniques on the results of spectral fits. We found that, 
apart from small residuals at the fluorescence emission line ener- 
gies, the background-subtracted spectra in each case were almost 
identical with the results of the spectral fits, being consistent to 
within the confidence range of the fit parameters. Hence, given 
that the blank-sky data more accurately accounted for the instru- 
mental background, we discuss the results obtained using this 
method. 



3.3.1. Spectral fits 

We extracted X-ray spectra from Region A and Region B (with 
point sources masked). Spectral anal ysis was performed using 
Xspec version 12.5.1pj( |Arnaud|l996l l. Initial fits on the regions 
showed that Region B is adequately fitted with either a soft 
collisional ionisation equilibrium model (CIE, the APEC model 
in Xspec, 



Smith et al. 



1 2001[ ), or a soft non-equilibrium ionisa- 
tion model (NEI, the NEI model in Xspec, Borkowski 2000), 
whereas Region A can be fitted with either a moderately hard 
single NEI model, a moderately hard NEI + soft CIE model, 
or a two NEI model, with the two-component models yielding 
the better fit statistics. The single-component model parameters 
of Region B were very similar to the softer component in 
the two-component fit to the Region A spectra. Therefore we 
conclude that the soft components in Regions A and B are due 
to the same diffuse source present in each region. Because of 
this and because of the faint, dispersed Ha structure surrounding 
Region B, we suspect that the soft X-ray emission in Regions A 
and B is due to a blowout region where hot gas is leaking from 



1 A detailed comparison of the different methods to estimate the 
backgrounds can be found in Sasaki et al. 1 20041. 

2 Note: Throughout the analysis spectra from both the MOS cameras 
for each observation were fitted simultaneously in Xspec. 



3 



P. J. Kavanagh et al.: XMM-Newton view of the N 206 superbubble 




30.0 32:00.0 30.0 31:00,0 30.0 3:30:00,0 30.0 29:00.0 32:00.0 30.0 31:00.0 30.0 9:30:00.0 30.0 29:00.0 32:00.0 30.0 31:00.0 30.0 5:30:00.0 30.0 29:00.0 

Right Ascension (J20D0) 



Fig. 3. (a) - MCELS Ha image with contours. The lower contour limit has been adjusted to highlight the Ha shell that confines 
the X-ray emitting gas. The levels correspond to four logarithmically spaced intervals between this adjusted lower limit and the 
maximum pixel value, (b) - 0.3-1 keV XMM-Newton image with the Ha contours overlaid. - (c) 0.3-1 keV XMM-Newton image 
with superbubble emission regions indicated. 



the superbubble interior into the surrounding regions. Given that 
the emission from the superbubble itself is best described by an 
NEI model, we favour this for describing the gas that escapes 
from the bubble and thus only the two component NEI model is 
considered for the remainder of the analysis. 

We refitted the spectra linking the softer NEI model parame- 
ters in Regions A and B to obtain better statistics on the soft com- 
ponent, while the harder NEI component was used to fit Region 
A only. The normalisations of the linked soft NEI components 
were adjusted to account for the difference in volume of the emit- 
ting regions. Because we assumed that the soft NEI emission 
comes from the hot gas in the blowout, the plasma density is 
uniform across both regions, only varying in emitting volume. 
Given that the model normalisation is defined as 



data and folded model 



1 



10 14 x 4nD 2 



J" n e n 



fgas dV, 



(1) 



where f gas is the filling parameter of the X-ray emitting gas, the 
ratio of the normalisations is that of the volumes. We assumed 
that the volume of the entire blowout (i.e., across Regions A 
and B) is defined by an ellipsoid 4' in height and 8' in width, 
determined from our X-ray images. Since only two-dimensional 
information is available for the blowout, we must assume its 
depth, which was set equal to its height of 4'. Thus, we have 
ellipsoidal radii of a = c — 2' and b — 4', or a — c w 28 pc and 
b « 56 pc at the distance of the LMC. To split the volume we cut 
the b axis of the ellipsoid into 5' and 3' divisions, corresponding 
to the blowout volumes in Regions A and B, respectively. This 
is not a perfect solution, but without additional information the 
true dimensions of the blowout are impossible to determine. 

The best-fit 2 NEI model is shown in Fig. [4] For clarity 
only the MOS 2 spectra from the effective 23 ks observation 
are shown. For completeness, the fit results for each of the 
models discussed are shown in Table [2] The absorbing h ydrogen 
column (A^h ) for the fits was fixed at 0.7 X 10 21 cm" 2 (Dickey 
|& Lockman| [T990) and the abundance was fixed at 0.5 solar. 
We chose to fix the A^h value for two reasons. Firstly, allowing 
the Nh parameter to vary results in a best-fit value lower than 
0.7 x 10 21 cm" 2 but consistent within the upper 90% confidence 




Energy {keV) 

Fig. 4. Best-fit 2 NEI model plotted on the MOS 2 spectra of 
the superbubble and blowout regions from the effective 23 ks 
observation. The spectrum of Region A is plotted in black, the 
spectrum of Region B in red. The components of the spectral 
fits are indicated by the dotted lines. The spectra are binned to a 
signal-to-noise ratio of 3. 



interval with the Dickey & Lockman ( 1990 1 value. Secondly, 
one could set the Na value to that of the background AGN 
[SHP2000] LMC 211 (see SectionO). However, with a derived 



value of ~ (1 + 2) x 10 21 c m" 7 , the Af H value for [SHP 2000] 
211, while consistent with the |Dickey & Lockman| ( fT9 90) value, 
is poorly constrained due to poor statistics. For these reasons we 
preferred to fix Nh at the known Galactic foreground value. 



4. Discussion 

4. 1 . Superbubble X-ray emission 

The morphological comparison between the X-ray and Ha 
images in Fig. [3] shows the brightest region of X-ray emission 
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CIE / NEI 






NEI 






Region (Model) 


(keV) 


(10 10 s cm- 3 ) 


norm i 
(10- 4 ) 


kT 2 
(keV) 


T2 

(10 10 s cm- 3 ) 


nor mi 

do- 4 ) 




A (NEI) 








1 01 140 


1 73 226 


1 09' 25 


219.11/171 = 1.28 


B (CIE) 
B (NEI) 


0.27 029 
79^ 

"•'^0.47 


9 106.03 

z " ly 0.99 


1.57'f 
25°^ 








46.20/56=0.83 
45.64/55=0.83 


A+B (CIE+NEI) 
A+B (NEI+NEI) 


25 27 

Q 24 

70 m 

/u 0.48 


1 74 3 - 64 


9 7^3.02 
z.. / j 2 43 

67° 59 


1 43 2 08 
99' 23 


1 52 219 
2.28 3 ;i 4 


67° 83 

1 29 1 59 


265.48/227=1.17 
262.18/226=1.16 



The upper and lower limits correspond to the 90% confidence intervals of the fit parameters. 



confined by a Ha shell structure, consistent with the superbubble 
picture of a hot diffuse gas confined by a cooler shell. However, 
in a rich stellar population such as that in the N 206 complex, 
an unresolved low-mass stellar population can account for a 
significant fraction of the seemingly diffuse emission. Indeed, 
Oskinova (2005) found that very early in the lifetime of a 
massive stellar population, the unresolved low-mass stellar 
population is the dominant contributor to the diffuse emission 
in the cluster. However, the stellar populations in the N 206 



stars for the remainder of the analysis. 



superbubble are at a stage of evolution (see Section 4.3.3 i where 
the truly diffuse emission is enhanced due to the higher mass 
loss rates of the post-main sequence population and SNRs. At 
this stage, the truly diffuse emission will dominate the apparent 
diffuse emission caused by the unresolved sources (Oskinova 
2005). To obtain an estimate for the contribution of unresolved 
low-mass population, we compared the N 206 stellar population 
to the well-known and characterised stellar population of the 
Orion Nebular Cluster (ONC), determined in the Chandra Orion 
Ultradeep Project (COUP, pitman et al.|2005| l, a method which 
has been used for a mu ltitude of stellar p opulation analyses 
(see |Getman et al.|[ 2006; Ez oe et al.|[2006l |Broos et al.||2007) 
|Kavanagh et al.|201 1| for example). We assumed that the N 206 
clusters and ONC have the same initial mass function (IMF), 
and X-ray luminosity distribution, differing only in the size 
of their underlying populations. This allowed us to scale the 
integrated X-ray luminosities of the lower mass bins of the ONC 
population to that of the N 206 population using the ONC mass 
stratified star counts and X-ray luminosities of Feigelson et al 



0.75 



(2005). Additionally, since the X-ray luminosity of youn g low- 
mass stars can decrease by as much as L x < 



( Preibisch 



& Feigelson 2005), the ONC X-ray luminosity distribution 
was corrected for the increased ages of the N 206 populations 
(determined in Section 4.3.3 I. We chose the lowest age limit for 



the N 206 populations for this correction, which corresponds to 
the period of highest X-ray luminosity for the low-mass stars. 
The high-mass stellar populations of the N 206 clusters was 
used to extrapolate the assumed IMF to lower masses. The 
X-ray luminosity of the ONC population was then scaled based 
on the number of stars in the 1-3 Mq mass bin. This provided 
us with an upper limit for the contribution of the unresolved 
low-mass population of < 4 x 10 34 erg s . Based on the spectral 
analysis in Section 3.3.1 we determined the X-ray luminosity 



of the N 206 superbubble to be ~ 7 x 10 J erg s . Hence, the 
unresolved low-mass population can at most contribute ~ 6 % to 
the observed X-ray emission. Considering that this upper limit 
is already quite low, we ignore the contribution of the low-mass 



As discussed in Section [372] no X-ray emission is detected 
from the majority of Region SB, which is unusual because 
we would expect the hot X-ray emitting gas to fill the entire 
superbubble volume. The likely explanation for this becomes 
apparent if we consider the spectra of the brightest point sources 
in Region A and in Region SB excluding Region A (which we 
will call Region SB -A), namely sources 4 and 6 identified as 
[SHP2000] LMC 211 and HD 269676, respectively. [SHP2000] 
LMC 211 is likely a background AGN, whereas HD 269676 
is a member of the NGC 2018 cluster in the superbubble. 
We performed a spectral analysis of these sources and found 
[SHP2000] LMC 211 to have a significantly lower absorbing 
hydrogen column (7Y H ~ (1 + 2) x 10 21 cm 2 ) than HD 
269676, whose corresponding value is several times higher 
(Nu ~ (9 ± 4) x 10 21 crrT 2 ). This indicates that the line-of-sight 
absorption to the Region SB -A is likely sufficiently high to 
absorb the soft superbubble X-ray emission. We cannot exclude 
the possibility that some of the line-of-sight absorption to HD 
269676 may be caused by stellar wind material in its circum- 
stellar environment and, thus, the Nh value of ~ (9 ± 4) x 10 21 
crrT 2 is an upper limit for the absorbing hydrogen column to 
Region SB-A. We note that for Galactic stars of similar spectral 
type this circumstellar abso rption can be as high as ~ 5 x 10 21 
crrT 2 (see Naze et al.|2011 for example). To estimate the effect 
that this additional absorption has on the superbubble X-ray 
emission from Region SB-A, we used the best-fit model of the 



detected X-ray emission from Region A (see Section 3.3.1 
added the additional absorption and considered the effect on 
the background-subtracted surface brightness of the emission. 
We find that if we assume no circumstellar absorption around 
HD 269676 (i.e., the superbubble X-ray emission is subject 
to the best-fit Nh), the surface brightness in the 0.3-1 keV 
range is expected to be almost an order of magnitude lower 
than the detected superbubble emission. If we consider that 
circumstellar absorption contributes ~ 5 x 10 21 crrT 2 , as is 
possible for Galactic stars of similar spectral type, then the 
surface brightness is still some 80% lower than that of the 
detected superbubble emission, and is lower than half of the 
blowout emission in Region B, which itself is quite faint. 



Additional evidence of significant variation in line-of-sight 
absorbing material comes from the Hi 21 cm line map towards 
N 206. We obtained the ATCA-Parkes 21 cm emission image in 
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the direction of N 206 from the H i Magellanic Clou d Surve>p] 
The ATCA and Parkes observations are outlined in lKim et al.l 
T998] > and |Kim et al] ( pOOl) , respectively. The Hi 21 cm line 



map (shown in Fig.|5J) clearly illustrates the variation in absorb- 
ing material across N 206. The Hi gas is morphologically quite 
similar to the distribution of the 24 /urn emitting dust in N 206 
( Gorjian et al.|2 004) with the same 'window' in the direction of 
Regions A and B. Additionally, the cooler dust, observed in the 
70 /im and 160 fim Spitzer images, shows the same variation in 
density across the superbubble. If we assume that the Hi gas in 
the N 206 region is distributed evenly between the foreground 
and background of the superbubble, the amount of absorbing 
material in Region SB-A is at least twice that of Region A itself. 
Indeed, estimation of the N H for [SHP2000] LMC 21 1 and HD 
269676 from the radio flux densities at their positions yields 
values of ~ 1 x 10 21 crrT 2 and ~ 4 x 10 21 ctrT 2 , respectively, 
which is roughly consistent with the best-fit A^h values to their 
X-ray spectra. The A^h value of HD 269676 is slightly lower 
than determined from its X-ray spectrum, but this is likely due 
to the assumption that the Hi gas its direction is distributed 
evenly between the foreground and background of the source. 
We note that not all of the H i gas in these considerations may 
be associated with the N 206 region. However, investigation of 
the velocity profile shows that the brightest and most significant 
21 cm emission is approximately centred around a heliocentric 
velocity of ~ 241 km s" 1 , which is the heliocentric velocity of 
the N 206 complex determined by Kim et aL] ( |2007) . Hence, we 
find it likely that the variation in absorption across the face of 
N 206 is sufficient to obscure the X-ray emission from much of 
the superbubble. 



4.2. Properties of the X-ray emitting gas 

The following calculation of the physical properties of the hot 
gas in the superbubble follows the method outlined in Sasaki 
|et al.| ( |201 1) . The superbubble emission is described by the 
harder of the two NEI components. From Equation [T] we can 
determine the gas density in the bubble. For the LMC metallic- 
ity £lmc = 0.5, n e = (1.2 + 0.013£lmc) gives n e w 1.20/1, with 
n being the hydrogen density. We must also consider the value 
of f gas , the filling parameter of the hot gas in the bubble and 
blowout. For a young stellar bubble, as is the case with N 206, 
the filling parameter is likely f gas 
[20041 



Sasaki et al. 201 1 



1 (seelOey & Garcia-Segura 
for example). H 



lence, we adopt this 

value of unity for the forthcoming calculations. Thus, Equation 
Q]can be written as 



1 



1.2n 2 V 



10 14 x4ttD 2 
and the gas density in the bubble is given by 



4.8 X 10 3l) x 



V 



cm 



-3, 



(2) 



(3) 



In the case of Region A, the volume is that of a cylinder-sphere 
intersection with the base radius of the cylinder equal to the 
radius of Region A, the radius of the sphere equal to the radius 
of Region SB, and the cylinder edge touching the sphere edge. 
We determine the volume of Region A (hereafter referred 
to as Va) to be ~ 8 x 10 60 cm 4 and, using the NEI model 
normalisation of norm^ = (1.3} g) x 10~ 4 , we estimate a gas 
density of n SB = (1.9 ± 0.2) x 10~ 2 cnr 3 . 



Using this gas density and the NEI model plasma tempera- 
ture kT2 = 0.99(±0.36) keV, the pressure of the gas inside the 
bubble, assuming the pressure is constant throughout the entire 
volume, is given as 



u -71:00:00.0 
o 




v A 



31:00.0 5:30:00.0 
RA 



Fig. 5. ATCA-Parkes 21 cm line map in the direction of N 206 
from the H i Magellanic Cloud Survey. The superbubble emis- 
sion regions (A and B) and superbubble boundaries (SB) are in- 
dicated. The contours correspond to 25%, 50%, and 75% of the 
peak emission. An additional region, shown here in blue, corre- 
sponds to the region from which we extracted the 21 cm line flux 



PsB/k = (n e + \.\n)T = 23n SB To = (5.1 ± 1.9) x 10 5 cm" 3 K, 

(4) 

With these derived values we can now estimate the thermal 
energy content inside the entire superbubble (Region SB), given 
by £sb(0 = 3/2P sb Vsb = (3.0 ± 1.1) X 10 51 erg, where V SB 
is the volume of the superbubble sphere with r = 4' (~ 56 
pc at the LMC distance), determined from our X-ray im- 
ages. In addition, the mass inside a homogeneous bubble is 
M SB = 2.31n SB XjUXra H xV SB = 599(±63)M , where fi = 0.61 
is the mean molecular weight of a fully ionised gas and mu is 
the hydrogen mass. 



to 



We applied a similar treatment to the 
gion using the softer NEI component in our 
our assumed blowout region volume 
density n BL = (1.3 ± 0.2) x 10~ 2 
P BL /k = (2.4 + 1.2) x 10 5 cnr 3 K, 
content E Bh (t) = (4.7 ± 2.4) x 10 50 erg., 
blowout region M B l = 142(+22) Mq. The overall mass and 
thermal energy content of the superbubble and blowout region 
are 741(+85) M Q and (3.5 ± 1.3) x 10 51 erg. 



blowout re- 
model with 
derive the gas 
cm" 3 , the pressure 
the thermal energy 
and the mass in the 



density for the analysis discussed in Section 4.3.2 



Available at http://www.atnf.csiro.au/research/HI/mc/ 



4.3. Superbubble energy budget 

The thermal energy stored in the X-ray emitting gas of the su- 
perbubble is supplied by the mechanical input of the stellar pop- 
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ulation. The pressure of the hot gas drives the expansion of the 
shell into the surrounding H i cloud. Logically we would expect 
that the thermal energy stored in the hot gas and the kinetic en- 
ergy of the shell should balance the mechanical energy input of 
the stellar population. However, observations of X-ray superbub- 
bles show this not to be the case. Rather, the observed energy 
stored in some superbubbles is much less than that injected by 
the stellar population (see |Cooper et al. 2004 |Maddox et al. 
2009 1 Jaskot et al. 201 1| for example). When observed super- 
bubble properties are compared to predictions from the standard 
pressure-driven model of Weaver et al. ( 1977[ ) or similar models, 
it has been found that the predicted dynamical evolution of the 
superbubble for a given stellar input proceeds at a slower rate 
than is observed. This is known as the growth rate discrepancy. 
It is clear that much of the energy injected into the superbub- 
bles by the stellar populations is 'missing' and several mech- 
anisms, from superbubble blowouts to non-thermal X-ray pro- 
duction, have been proposed to explain this problem. Detailed 
discussions on this subject can be found in Cooper et al. (2004) 
and Jaskot et al. (2011). We now examine N 206 with regard 
to the growth rate discrepancy, determining the kinetic energy 
stored in the shell and, as the ionisation front may be trapped in 
the shell, the surrounding H i gas, as well as the energy supplied 
by the stellar population. 



should be considered an upper limit as not all of the H i gas in 
the region may be associated with the superbubble. 

4.3.3. Stellar content and input 

To date, the high-mass stellar population of the N 206 Hn 
region has not been the subject of a detailed spectroscopic 
analysis. However, the photometric study of high-mass stars in 
the Magellanic Clouds of Massey et al. ( 1995| ) provides us with 
the spectral types of the massive stars in N 206. Additionally, 
more recent analyses of specific source classes in the LMC 
(Breysach er et al.| |T999; Rob ert et al.||2 003) supplement and/or 
update the spectral types of [Massey et al. (1995). We used the 
SIMBAD database to determine the most recent information 
available for the OB stars in the N 206 region and adopted 
the determined spectral types. Association membership was 
assigned based on the location of the stars with respect to the 
association dimensions for LH 66, LH 69, and NGC 2018, as 
determined by Bic a et aL[ ([1999 ). Since we only have limited 
data on these sources, it is difficult to say whether they indeed 
belong to the stellar associations or just appear to be due to 
projection effects. Nevertheless, we list the high-mass stars with 
known spectral types with their likely host association in Table[3] 



To determine the mass and energy input of the massive stars, 
we followed a similar method to that of Sas aki et"aX] ( |201 1). We 
assigned effective temperatures and luminosities from Martins 
[eTaLl ( |2005] ) for the O stars and |Schmidt-Kaler| ( fT982> for the B 
stars, and used the mass-luminosity relation of Vitrichenko et al. 
(2007) to obtain the mass (M) of each star: 



19(M/M q ) 2J6 Lq. 



4.3.1 . Energy stored in the Ha shell 

To determine the kinetic energy stored in the Ha shell, we must 
first determine its mass. This can be derived from the Ha lu- 
minosity of the shell, L Ha . Using scans of the Curtis Schmidt 
plates of |KennicutT & Hodge ( fT9§6| ), |Dunne et al.| ( |200l) deter- 
mined the shell Lu a to be ~ 9 x 10 erg s" 1 . For ionised gas at 
temperatures of ~ 10 4 K, L Ha = 3.56 X lO' 25 n e n p Vf she u, where The radii of the stars are obtained from 
n e and n p are the number densities of electrons and protons, 
respectively, V is the volume and f s i ie n is the filling parameter. 
Helium is likely singly ionised, so n e = \.\n p . Substituting this 
into the equation above with the known value of V, we obtain 

— 1/2 ^ 

n e — \.Q2f sl J n cm . The total mass in the superbubble shell is 

then given by M sh = l.21n e Vf she ,m H * 3 x 10 4 M f 1 ^. The 
filling parameter fsheg is determined from the thickness of the 
Ha shell, which in the case of our superbubble is ~ 10 - 15% of 
the superbubble radius. This corresponds to a filling parameter 
f shell = Q.35(+0 05), yielding a tota l mass of (1.8+0.2)xl0 4 M Q . 
We know from IDunne et al.l (120011) that the expansion velocity 



R 



(5) 



(6) 



where T e g is the effective temperature of the star and cr is the 
Stefan-Boltzmann constant. The wind velocity was then calcu- 
lated using the theory of radiation-driven winds ( |Castor et al. 
[1975]), 



2GM 



R 



X (1 - L/Ledd) 



0.5 



(7) 



of the Ha shell is ~ 30 km s , thus we derive a kinetic energy 
E k = \Mv 2 * (1.6 ± 0.2) x 10 50 erg. 



4.3.2. Energy stored in H i gas 

The mass and energy input of the stellar populations in the super- 
bubble drive the shell into the surrounding Hi gas. If the ionisa- 
tion front of the stellar population becomes trapped in the shell, 
the outer layer remains cold Hi gas. The mass of this swept-up 
H i gas can be determined from the 21 cm line emission line data 
from the ATCA-Parkes survey of the Magellanic Clouds, used 
to produce Fig. [5] We determine the total line flux density from 
the region of the superbubble and blowout (shown as the blue 
region in Fig. [5]) to be ~ 238 Jy. This corresponds to a mass of 
~ 1 x 10 5 M at the LMC distance of 48 kpc. We assumed that 
the expansion velocity of the swept up H i gas is the same as that 
of the Ha shell, namely, ~ 30 km s _1 . This implies the kinetic 
energy stored in the H i gas is ~ 1 x 10 51 erg. However, this value 



where a = 2.6(+0.2) ( |Lamers et al. 1995) ), v esc is the photo- 
spheric escape velocity and L e dd is the Eddington luminosity 
given by L e[ jd = 4-nG xMxm p xc/<TT with m p being the proton 
mass and cr T the Thompson cross section for the electron. The 
mass loss is determined from the single scattering limit 



M : 



x c 



(8) 



We corrected the wind velocities and mass loss rates for LMC 
metallicity following Leith erer et al.| ( | 1992 ) and Vink et al. 
(2001), respectively. The estimates for wind velocities and mass 
loss rates, which are shown in Table [3] are in general agreement 
with the wind properties of LMC stars in the literature (see 
Moki errTet al.|2007| for example). 

To determine the total mass and energy injected into the 
superbubble by the stellar populations we must first determine 
their ages. However, because no high-quality optical data are 
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Table 3. High-mass stellar content of superbubble. Voo and M are the estimated terminal wind velocity and mass loss rates, respec- 
tively (see text) 



RA 


DEC 


Object Name 


Association 


Spectral Type 


Voo 


M 


















(kms- 1 ) 


(Mq yr" 


) 


05 


29:52 


-71 


04:07 


[MLD95] LMC 1-548 


LH 66 


B1.5V 


1790 


5 A ss 
JM X 


1 u 


-8 


05 


29:52 


-71 


04:33 


2MASS J05295152-7104325 


LH 66 


08V 


2094 


A 1. v 
■4. J X 


1 u 


-7 


05 


30:11 


-71 


04:04 


2MASS J05301 149-7104034 


LH 66 


B1III 


1598 


Z.o X 


i fr 


-7 


05 


30:11 


-71 


04:47 


2MASS J05301 122-7104464 


LH 66 


09V 


2039 


9 Q v 
L.y X 


1 fl" 

i u 


-7 


05 


30:12 


-71 


06:00 


2MASS J05301 167-7105596 


LH 66 


B0.5V 


1872 


1 Q v 

L.y X 


1 u 


-7 


05 


30:19 


-71 


06:07 


2MASS J05301948-7106071 


LH 66 


B0V 


1942 


^ 1 V 
j . L X 


i fl - 
1 u 


-7 


05 


30:23 


-71 


06:19 


2MASS J05302268-7106194 


LH 66 


08V 


2066 


A A v 


i n- 
1 u 


-7 


05 


30:25 


-71 


03:45 


2MASS J05302469-7 103445 


LH 69 


07V 


2081 


7 a v 
/ .U X 


1 fl" 
1 u 


-7 


05 


30:31 


-71 


02:32 


2MASS J05303066-7102315 


LH 69 


B0.5V 


1872 


1 Q ss 
L.y X 


1 fl" 
1 u 


-7 


05 


30:35 


-71 


01:57 


2MASS J05303520-7101570 


LH 69 


BlVe 


1829 


1 fl v 
1 .U X 


1 fl" 
1 u 


-7 


05 


30:37 


-71 


01:43 


[MLD95] LMC 1-715 


LH 69 


B0.5V 


1872 


1 Q v 
L.y X 


i fl - 
1 u 


-7 


05 


30:37 


-71 


04:16 


2MASS J05303660-7104158 


LH 69 


09V 


2039 


x 


i fl - 
1 u 


-7 


05 


30:39 


-71 


01:48 


HD 37248 


LH 69 


WC4+09 


2400 


J.J X 


i fl - 
1 u 


-5 


05 


30:45 


-71 


01:45 


[MLD95] LMC 1-713 


LH 69 


B1III 


1598 


Z.o X 


1 fl" 
1 u 


-7 


05 


30:47 


-71 


03:16 


2MASS J05304707-7103156 


LH 69 


BOIa 


1369 


9 9 y 


1 0" 
1 u 


-6 


05 


30:48 


-71 


04:02 


HD 269660 


LH 69 


Bllab 


1189 


1 S V 

1 .J X 


1 fl" 
1 u 


-6 


05 


30:51 


-71 


02:48 


2MASS J05305063-7 102480 


LH 69 


08V 


2034 


S 9 v 
J.Z, X 


1 fl" 
1 u 


-7 


05 


30:54 


-71 


05:11 


2MASS J05305448-7105105 


LH 69 


Bill 


1598 


2.8 x 


10" 


-7 


05 


30:58 


-71 


01:32 


2MASS J05305766-7101315 


LH 69 


B1III 


1598 


2.8 x 


10" 


-7 


05 


30:58 


-71 


03:34 


2MASS J05305760-7 103337 


LH69 


08V... 


2066 


4.4 x 


10" 


-7 


05 


31:10 


-71 


03:52 


2MASS J05310966-7103515 


LH69 


B1III 


1598 


2.8 x 


10" 


-7 


05 


31:11 


-71 


01:43 


2MASS J05311070-7101433 


LH69 


08III... 


1788 


1.3 x 


10" 


-6 


05 


31:12 


-71 


02:58 


[MLD95] LMC 1-485 


LH69 


BOVe 


1942 


3.1 x 


10" 


-7 


05 


31:12 


-71 


04:10 


[MLD95] LMC 1-550 


NGC 2018 


05V 


2095 


1.8 x 


10" 


-6 


05 


31:14 


-71 


03:51 


[MLD95] LMC 1-576 


NGC 2018 


07.5V 


2074 


5.6 x 


10" 


-7 


05 


31:16 


-71 


04:10 


HD 269676 


NGC 2018 


04-5III(f) 


1802 


3.2 x 


10" 


-6 


05 


31:16 


-71 


03:43 


[MLD95] LMC 1-552 


NGC 2018 


04Iab:... 


1435 


8.0 x 


10" 


-6 


05 


31:19 


-71 


04:39 


2MASS J05311887-7104389 


NGC 2018 


07Ve 


2081 


7.0 x 


10" 


-7 


05 


31:22 


-71 


04:06 


2MASS J05312153-7104057 


NGC 2018 


08V 


2066 


4.4 x 


10" 


-7 


05 


31:24 


-71 


04:13 


[MLD95] LMC 1-546 


NGC 2018 


06Ve 


2094 


4.2 x 


10" 


-6 



available, robust age determination methods such as theoretical 
isochrone fitting are beyond our grasp. One possibility is to 
assume the age of the superbubble as the age of the stellar 
populations. |Dunne et al. ( 200 1 [ ) gave the expansion velocity 
of the superbubble Ha shell as ~ 30 km s _1 . Because the 
superbubble has expanded to a diameter of ~ 112 pc, this 
suggests an age of ~ 2 Myr. However, the expansion veloci ty of 
~ 30 km s falls in the high-velocity superbubble regime (Oey 



and luminosities at Galactic metallicit)^] we used the theoretical 



1996 ), meaning the N 206 superbubble must have undergone 
induced acceleration at some point in its history, most likely 
due to SNRs impacting the outer shell. Thus, the superbubble 
is likely older than 2 Myr and, as such, is not a good estimate 
for the age of the stellar populations. Gorjian et al. ( 2004| l used 
previous radio and optical photometric measurements of |Kim| 
et al. ( 1999 1 and |Bica et al.| ([T996), respectively, to constrain the 
age of the stellar populations to 2-10 Myr. However, this wide 
range of possible ages does not allow the total mass and energy 
input of the stellar population to be sufficiently constrained. 
Instead we took the effective temperatures and luminosities used 
above to produce HR diagrams for each association to compare 
them with theoretical isochrones. While this is not a perfect 
solution, it at least allows us to obtain better age estimates than 
are currently available. Since we assigned effective temperatures 



isochrones of Schaller et al. (1992| at z — 0.020 to gauge the 



likely ages of the stellar associations. The HR diagram is shown 
in Fig. 

The age of the NGC 2018 population is well described 
by the 2 Myr isochrone. Owing to slight scatter in the points 
and the magnitude of the error bars, we assign an age range of 
2(±1) Myr for NGC 2018. With regard to LH 66 and LH 69, 
the presence of the Bl III stars, which indicate an age > 10 
Myr, is somewhat at odds with the main-sequence populations, 
which are indicative of much younger ages. Because of this, 
we suggest that the Bl III sources are not associated with 
LH 66 or LH 69. This is a reasonable conclusion because the 
wider N 206 Hn region hosts stellar populations outside of 
the superbubble we discuss in this paper. Due to dynamical 
ejection or supernova kick processes, a star with a modest 
ejection/kick velocity of ~ 30 km s _1 can cross the diameter 
of the superbubble (~ 112 pc) in < 5 Myr, well within the 
lifetime of early B-type stars. The remaining LH 69 stars show 
a reasonably large spread in age but the location of the 08 III 
and B0 la stars suggest an age of ~ 5 Myr. The Bl lab star 
is older than this while the top of the main sequence suggests 



4 We corrected the results of our mass loss and terminal velocity cal- 
culations for metallicity effects and not the input effective temperatures 
and luminosities, which are representative of solar metallicity. 
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Fig. 6. HR diagram showing the stellar populations of LH 66, LH 69, and NGC 2018 compared to the evolutionary tracks and 
isochrones of Schaller et al. ( 1992 ). LH 66 sources are indicated by the blue inverted triangles, LH 69 sources by the green triangles, 
and NGC 2018 sources by the red circles. The solid black line is the ZAMS isochrone with the dashed lines showing isochrones 
from 1-10 Myr in increments of 1 Myr. The dotted lines show the evolutionary tracks of stars with masses of (from top to bottom) 
120 Mq, 85 Mq, 60 Mq, 40 M©, 25 Mq, 20 M©, 15 M©, and 12 M©. For clarity, we only show part of the evolutionary tracks of 
the most massive stars. We also omit the main-sequence B stars because they are unimportant in this age determination. The spectral 
types of the post-main sequence stars are also labelled to aid in the discussion. Uncertainties are unavailable for the B stars. 



an age younger than 5 Myr. Therefore, we set a cautious age 
range of 5(+3) Myr. Unfortunately, the LH 66 population has 
no post-main sequence objects. We assigned its age based on its 
upper main sequence to be 4(+2) Myr. The age determinations 
for each of the stellar associations are consistent with the 2-10 
Myr age constraint provided by Gorjia net al.| ( [20 04). However, 
our age estimations consider only single-star evolution and not 
close binary evolution effects. Vanbeveren et al. ( 1998 ) showed 
that the majority of massive stars are formed in interacting close 
binary systems and that mass transfer effects in these binaries 
can lead to incorrect age determinations for massive stellar 
clusters older than ~ 4 Myr through to the process of 'starburst 
rejuvenation'. This leads to a possible underestimation of 
cluster age of up to a few Myr if close binary evolution is not 
considered. For the N 206 associations with ages > 4 Myr, 
namely LH 66 and LH 69, the determined age ranges are large 
due to the limitations of the available data, but may sufficiently 
account for this underestimation of cluster age. This question 
can only be resolved with higher quality optical/IR data and the 
application of binary evolutionary tracks. 



Using these determined ages, mass loss rates, and wind 
velocities, we estimate that the O and early-B stars in the 
superbubble have supplied a mass of 86(±52) M© and a total 
energy of (2.5 + 1.5) X 10 51 erg. 

The single Wolf-Rayet (WR) star in the superbubble 
HD 37248, which is of carbon class (WC) and a member of LH 
69, also supplies a large amount of mass and energy. According 
to the theoretical evolutionary tracks at LMC metallicity of 
Schaerer et al. ( 1993| l, only a star more massive than 60 M© will 
enter the WC phase. Given the large error on the age estimate 
for LH 69, namely 5(±3) Myr, it is impossible to constrain 
the initial mass of the WC star further. We adopted mass loss 
rates and velocities from Leither er et al.| ([1997 ) and corrected 
the mass loss rates for the metallicity of the LMC according 
to |Crowther| ( |2007| ) (listed in Table [3]), the results of which 
are in rough agreement with those of other WC class stars in 
the LMC from |Grafener efaLl ([1998). We calculated the total 
mass and energy input over the WR lifetimes (from Schaerer 
et al.| 1993] > for stars of mass > 60 Mq, assuming they are half 
way through their WC lifetime. Additionally, we calculated 
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the mass and energy input of the star before it entered the WR 
phase, which we estimated using the methods outlined above 
and the pre-WR lifetimes of stars > 60 M0 from| Schaerer et aL] 
(1993). We find that over its lifetime the WC star has injected 
20(+3)M o and (1.1 ± 0.2) x 10 5I erg into the superbubble. The 
mass loss rates and wind velocity estimates are determined for 
the evolution of a single WC star. Because HD 37248 is in a 
binary system, its age estimate and wind parameters are more 
uncertain. For example, Sander et al. ( 2012| i showed that stellar 
wind parameters of binary WC stars are different from single 
WC stars of the same spectral type (see their Fig. 6). However, 
without high-quality optical spectra we cannot further constrain 
the wind parameters of this object. 

We must also estimate the contribution of previous supernovae 
(SNe). We assumed that the stellar populations formed from 
the same parental molecular cloud and are described by a 
standard Salpeter IMF ( |Salpeter|1955] >. We find from Fig.[6]that 
there are some 16 stars in the 20-40 Mq bins and 4 stars with 
masses > 40 M0, including the WR star. Hence, we determined 
that 3(+2) SN have already occurred in the bubble. Taking an 
SN explosion energy of 10 51 erg, these SNae have injected 
~ (3 + 2) x 10 51 erg into the superbubble and, assuming an 
average initial mass of 60 M0, supplied some 60 - 300 M0 of 
material through their ejecta. 

Thus, in total the stellar population has injected 286(±175)M0 
and (6.6 + 3.7) x 10 51 erg into the bubble. 



4.3.4. Observed/expected energy and mass comparison 

The observed and expected energies and masses in the N 206 
superbubble are summarised in Table |4] In total, the overall 
energy input by the stellar populations is (6.6 + 3.7) x 10 51 erg 
while the observed energy stored in the superbubble is 
< (4.7 ± 1.3) X 10 51 erg. Hence, no definitive conclusions can be 
drawn regarding the energy budget in the N 206 superbubble. 
Because of the limited data on the stellar populations, the 
uncertainty on the age estimates of the stellar associations 
is substantial, making the uncertainty on the stellar energy 
input quite large. Only when high-quality optical/NIR data are 
available for N 206 will the stellar energy input be sufficiently 
constrained to say for sure. 

The overall mass contained in the superbubble and blowout 
region is 741 (+85) Mq. We expect the stellar populations to 
have contributed 286(±175)M0. Even with the large uncertain- 
ties on the mass from the stellar populations arising from the 
large uncertainty in ages, the observed mass is much higher. This 
higher-than-expected observed mass is likely due to mass load- 
ing by evaporation of entrained interstellar clouds and/or turbu 
lent mixing of material from the cold shell ( Weaver et al.|1977" 
|Silichetal.|1996l l. 



4.3.5. Evidence for an old SNR in N 206? 

During our analysis of the N 206 superbubble, we attempted 
to identify signatures of SNR impacts of the superbubble shell 
walls using the well-known optical identification method of the 
[S n]/Ha ratio, with a value of [S n]/Ha > 0.4 being indicative of 
a shock-heated nebula (Mathewson & Clarke 1973a; [Raymond 
|1979| l. However, as found by |Fesen et al.| ( 1985| and refer- 



Table 4. Observed/expected energies and masses in the N 206 
superbubble. Observed values are calculated assuming uniform 
gas density and a filling factor of unity. 





Energy 


Mass 




(10 51 erg s- 1 ) 


(Mq) 


Observed 


X-ray emitting gas 


3.5(±1.3) 


741(±85) 


Kinetic energy of Ho- shell 


0.16(+0.02) 




Kinetic energy of surrounding H 1 gas 


~ 1 




Total observed 


4.7(±1.3) 


741(±85) 


Expected Input 


Stellar wind input 


3.6(±1.7) 


106(±55) 


Supernovae input 


3(±2) 


180(±120) 


Total input 


6.6(±3.7) 


286(±175) 



ences therein), if an SNR is located in an Hn region, a higher 
[S n]/Ha > 0.67 is required to unambiguously attribute the op- 
tical emission to an SNR, particularly in the absence of addi- 
tional SNR identification criteria such as X-ray or radio emis- 
sion. To assess the [S n]/Ha ratio in the N 206 region, we used 
the MCELS data to create continuum-subtracted [S n] and Ha 
images, and produced an [S n]/Ha ratio image. Contours indi- 
cating the regions of N 206 with [S n]/Ha > 0.67 were de- 
rived from this image and are plotted on the MCELS Ha in 
Fig. [7] We found no evidence for SNRs in the X-ray super- 
bubble. Interestingly, however, we find a much larger, approx- 
imately semi-circular shell of shock heated material in the wider 
N 206 complex (highlighted as the dashed green circle in Fig.[7]i, 
as well as SNR B0532-71.0 in the east and the candidate SNR 



[HP99] 1234 (Haberl & Pietsch 1999 ) in the west. A study of the 
candidate SNR [HP99] 1234 will be published in a separate pa- 
per (fKavanagh et al.|2012[ ). Intriguingly, the known HMXB and 
candidate HMXB in the complex are both located close to the ap- 
parent centre of the large optical shell (see Fig. [7]). We determine 
the diameter of the shell to be ~ 195 pc. This is certainly large 
for an SNR, considering the largest confirmed SNRs in the LMC, 
such as SNR 0450-70.9 dMafhewson et al.|1985l|Williams et al. 



2004 



1976 



Cajko et al. 2009 ), L MC SNR J0550-682 3lDavies et al 



Filipovic et a 



6542(Klimek et al. 



. 1998| |Bozzettoet al.|2012 
2010)1, are significantly sma 



diameters can be ac 



and SNR0506- 
ler. Larger SNR 



:hieved if the ejecta expand into a low-density 
medium as is the case for the galactic Monogem Ring. Plucinsky 
|et al.| ( |1996| ) determined the diameter of the Monogem Ring to 
be ~ 133 pc and showed that this SNR has evolved in a region 
of the Galaxy with abnormally low density. These two exam- 
ples point to a possible explanation for the large diameter of 
the optical shell. We suggest that, in the case of the large op- 
tical shell, the ejecta have been expanding into a wind-blown 
bubble, likely a superbubble, and are only now impacting and 
shock-heating the superbubble shell walls, producing the [S 11] 



enhancement we observe. As discussed by Chu (1997i, this is 
an entirely plausible scenario. An SNR deep in a superbubble 
will not be detectable until the ejecta shock dense material, such 
as the shell walls or cloudlets in the superbubble interior. It is 
not unreasonable to assume the presence of a second older su- 
perbubble in the N 206 complex. As mentioned in our analysis, 
the wider complex hosts other OB associations. The association 
BSDL1927 ( |Bica et al.||1999| ) is located outside of the super- 
bubble discussed above, but inside the large optical shell. Using 
the SIMBAD database we determine that this association con- 
tains a small number of later O stars and a WR star. Indeed, 
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-70d50m00.0s 



71d00m00.0s 



I 
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Right Ascension (J2000) 

Fig. 7. MCELS Ha image of the N 206 H n region. Left: N206 superbubble and blowout regions shown with the white contours 
indicating regions with [Sn]/Ha > 0.67. Right: Known HMXB XMMU J052952.9-705850 and candidate HMXB USNO-B1.0 
0189-00079195 indicated by the small circles labelled 1 and 2, respectively. The large dashed green circle marks the apparent shell 
of the candidate old SNR. Again the white contours indicate regions with [S n]/Ha > 0.67. Both SNR B0532-71.0 and the candidate 
SNR [HP99] 1234 are evident in the east and west, respectively. 



the known HMXB, XMMU J052952.9-705850, is located ~ 2' 
east of this association. Unfortunately, the XMM-Newton obser- 
vations do not cover the west of the large optical shell where 
the shocks are located (see Pig. [7j so it is uncertain whether the 
shell is X-ray bright. Only high-quality X-ray and radio data will 
resolve the physical nature of this object. 



5. Summary 

Using XMM-Newton observations, we have detected significant 
soft diffuse X-ray emission in the N 206 Hn region. We 
performed a morphological comparison of the X-ray emission 
to Ha emission in N 206 utilising MCELS data and found that 
the brightest X-ray emission is confined by an Ha structure that 
extends as an approximately circular shell. This points to a su- 
perbubble origin for the detected emission where a hot shocked 
X-ray emitting gas is confined by the cooler swept-up material, 
confirming the results of |Dunne et al.| ( [2b01 1. In addition, faint 
diffuse X-ray emission was detected extending beyond the Ha 
shell, which we attribute to a blowout region where hot gas 
from the superbubble is escaping into the surrounding ISM. We 
performed spectral analyses on the superbubble and blowout 
emission, finding both to be thermal in origin, consistent with 
a hot shocked gas. We also used the MCELS data to identify 
a large semi-circular shell of shock-heated material with the 
optical emission line properties of an SNR. Because of the 
large diameter (~ 195 pc), we suggest it is likely an old SNR 
that has expanded into a wind-blown bubble and is only now 
interacting with the bubble walls to produce the optical emission 
we observe. 



Using the results of the N 206 superbubble spectral 
analysis and assumed volumes, we calculated the physi- 



cal properties of the hot gas in the bubble and blowout. 
We determined a temperature of kT = 0.99(+0.36) keV, 
a density of «sb = (1.9 ± 0.2) x 10 -2 crrT 3 and a pres- 
sure of P SB /k = (5.1 ± 1.9) x 10 5 cm" 3 K for the 
hot gas in the superbubble and corresponding values of 
kT = 0.70(±0.31) keV, n BL = (1.3 ± 0.2) x 10~ 2 cnr 3 and 
PBh/k = (2.4 + 1.2) x 10 5 cirT 3 K for the blowout. Using these 
values, we determined the overall mass and thermal energy 
content of the superbubble and blowout region to be 741(±85) 
M Q and (3.5+ 1.3) x 10 51 erg. 

To determine if the N 206 superbubble exhibits the same 
growth rate discrepancy as observed in other superbubbles, we 
compared the observed energy stored in the superbubble to the 
expected energy supplied by its stellar population. We combined 
the thermal energy stored in the hot gas with the kinetic energy 
of the expanding Ha shell and surrounding H i gas to determine 
a total observed energy of < (4.7 + 1.3) x 10 51 erg. We used the 
available information of the N 206 stellar population to deter- 
mine the expected energy from stellar winds and past SNRs of 
(6.6 ± 3.7) x 10 51 erg. Due to the poorly constrained stellar pop- 
ulation input we cannot definitively confirm a growth rate dis- 
crepancy for the N 206 superbubble. Only detailed optical/NIR 
analysis of the stellar population in the superbubble will con- 
strain the input sufficiently to provide a conclusive answer. 
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